Josephson effects in dilute Bose-Einstein condensates 
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We propose an experiment that would demonstrate the 
"dc" and "ac" Josephson effects in two weakly linked Bose- 
Einstein condensates. We consider a time-dependent barrier, 
moving adiabatically across the trapping potential. The phase 
dynamics are governed by a "driven-pendulum" equation, as 
in current-driven superconducting Josephson junctions. At 
a critical velocity of the barrier (proportional to the critical 
tunneling current), there is a sharp transition between the 
"dc" and "ac" regimes. The signature is a sudden jump of a 
large fraction of the relative condensate population. Analyti- 
cal predictions are compared with a full numerical solution of 
the time dependent Gross-Pitaevskii equation, in an experi- 
mentally realistic situation. 
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The Josephson effects (JE's) are a paradigm of the 
phase coherence manifestation in a macroscopic quan- 
tum system [0-|3|. Observed early on in superconduc- 
tors ||, JE's have been demonstrated in two weakly 
linked superfiuid 3 He-B reservoirs J5J. Weakly interacting 
Bose-Einstein condensate (BEC) gases || provide a fur- 
ther (and different) context for JE's. Indeed, magnetic 
and optical traps can be tailored and biased (by time- 
dependent external probes) with high accuracy al- 
lowing the investigation of dynamical regimes that might 
not be accessible with other superconducting/superfluid 
systems. The macroscopic BEC's coherence has been 
demonstrated by interference experiments [^,0, and the 
first evidence of coherent tunneling in an atomic array, 
related to the "ac" JE, has been recently reported ||. 

A superconducting Josephson junction (SJJ) is usu- 
ally biased by an external circuit that typically includes 
a current drive I ex t- The striking signatures of the 
Josephson effects in SJJ are contained in the voltage- 
current characteristic (V-I ex t), where usually one can 
distinguish between the superconductive branch or "dc"- 
branch (V = 0, I ex t 7^ 0), and the resistive branch or 
"ac" -branch (V « RI ex t) (see for example [||). Exter- 
nal circuits and current sources are absent in two weakly 
linked Bose condensates and the Josephson effects have 
been related, so far, with coherent density oscillations be- 
tween condensates in two traps or between condensates 
in two different hyperfine levels This collective 

dynamical behavior is described by a non-rigid pendu- 
lum equation g], predicting a new class of phenomena 
not observable with SJJ's. 

Now the following question arises: can two weakly 



linked condensates exhibit the analog of the voltage- 
current characteristic in SJJ? Although BECs are obvi- 
ously neutral, the answer is positive. A dc current-biased 
SJJ can be simulated by considering a tunneling barrier 
moving with constant velocity across the trap. At a crit- 
ical velocity of the barrier a sharp transition between 
the "dc" and "ac" (boson) Josephson regimes occurs. 
This transition is associated with a macroscopic jump 
in the population difference, that can be easily moni- 
tored experimentally by destructive or non-destructive 
techniques. 

In the following we will briefly introduce the phe- 
nomenological equations of the resistively shunted junc- 
tion (RSJ) model for the SJJ. We will describe the cor- 
responding experiment for two weakly linked BECs and 
show that the relevant equations are formally equivalent 
to the RSJ equations. Then we compare the analytical re- 
sults with a numerical integration of the Gross-Pitaevskii 
equation in a realistic 3D setup. 

In the RSJ model, SJJ is described by an equivalent 
circuit Q in which the current balance equation is 

I c sin{6) + GV + CV = I ext (1) 

where I c is the upper bound of the Josephson supercur- 
rent I (which is represented, in the ideal case, by the 
sinusoidal current-phase relation / = J c sin(#)); G is an 
effective conductance (offered by the quasiparticles and 
the circuit shunt resistor), and C is the junction capac- 
itance. The voltage difference V across the junction is 
related to the relative phase 9 by 

6 = 2eV/h . (2) 

In the low conductance limit G <C oj p C where lu p — 
y // 2eI c /hC is the Josephson plasma frequency, combining 
equations (|]J) and (@) leads to the "driven pendulum" 
equation 

e = -^U(9) (3) 

where U is the tilted "washboard" potential: 

U (0) = 1 - cos(0) + i6 (4) 

with i = Iext/ 1 c- This equation describes the tran- 
sient behavior before the stationary dissipative behavior 
is reached (resistive branch). If we start from equilib- 
rium, with i = 0, and increase adiabatically the current, 
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no voltage drop develops until the critical value i = 1 is 
reached (neglecting secondary quantum effects). At this 
point V continuously develops until a stationary asymp- 
totic dissipative behavior is reached in a time scale ap- 
proximately of order C/G. Similar phenomenology may 
occur in BECs and we will derive equations formally iden- 
tical to Equations (||) and (||). 

A weak link between two condensates can be created 
by focusing a blue-detuned far-off-resonant laser sheet 
into the center of the magnetic trap The weak link 
can be tailored by tuning the width and/or the height of 
the laser sheet. Raman transitions between two conden- 
sates in different hyperfine levels provide a different weak 
link [Q, in analogy with the "internal Josephson effect" 
observed in 70s with 3 He — A JTsJ . 

Here we consider a double well potential in which the 
laser sheet slowly moves across the magnetic trap with 
velocity v (but our framework can be easily adapted to 
investigate the internal Josephson effect). In the limit of 
very low v, the two condensates remain in equilibrium, 
i.e. in their instantaneous ground state, because of the 
non-zero tunneling current that can be supported by the 
barrier. In fact, an average net current, proportional to 
the velocity of the laser sheet, flows through the barrier, 
sustained by a constant relative phase between the two 
condensates. This keeps the chemical potential difference 
between the two subsystems locked to zero, as in the 
SJJ dc-branch. However, the superfluid component of 
the current flowing through the barrier is bounded by a 
critical value I c . As a consequence there exists a critical 
velocity v c , above which a non-zero chemical potential 
difference develops across the junction. This regime is 
characterized by a running-phase mode, and provides the 
analog of the ac-branch in SJJ's. 

The "dc" and "ac" BEC regimes are governed by a 
phase-equation similar to the current-driven pendulum 
equations ((H) and (^). Such equations together with the 
sinusoidal current-phase relation I = I c sin(f?) describe 
the phase difference and current dynamics. The dimen- 
sionless current i is related to the barrier velocity by 
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[H (t)+g\*\ 2 ]* 



(7) 



i = v I v c 

with the critical velocity v c given by 

= frj 

Vc ~ F 



(5) 



(6) 



where F is to a good approximation represented by dou- 
ble the average force exerted by the magnetic trap on 
single atoms in one well. 

Equations (^|)-(^|) can be derived by a time-dependent 
variational approximation and have also been verified, as 
we discuss below, by the full numerical integration |L8[ of 
the Gross-Pitaevskii equation (GPE) The GPE 

describes the collective dynamics of a dilute Bose gas at 
zero temperature: 



where Hq (t) = — ^ V 2 + V ex t (r, t) is the non interacting 
Hamiltonian and where g = 4:Trh 2 a/m, with a the scat- 
tering length and m the atomic mass. The order param- 
eter "J = \V (r,t) is normalized as / dr \\V (r,i) | 2 = N, 
with N the total number of atoms. The external po- 
tential is given by the magnetic trap and the laser bar- 
rier V ext {r,t) = V t ra P (r) + Vi aser (z,t). We consider 
a harmonic, cylindrically symmetric trap Vt rap (r) = 
\mtjj 1 (x 2 + y 2 ) + ^mui 2 z 2 where u) r and u>o are the radial 
and longitudinal frequency, respectively. The barrier is 
provided by a Gaussian shaped laser sheet, focused near 
the center of the trap Vi aS er (z) — Vo exp (— (z — l z ) 2 /X 2 ) 
with the coordinate l z (t) describing the laser motion and 
v = dl z /dt its velocity. 

The equations (||) to @ can be derived by solv- 
ing variationally the GPE using the ansatz: (r, t) — 
ci(t) ipi {r) + c 2 (t)ip 2 (r), where c„ = \/N n (t) exp {iO n (t)) 
are complex timc-dcpcndcnt amplitudes of the left n = 1 
and right n = 2 condensates (see also || ) . The trial wave 
functions ty\ t i (r) are orthonormal and can be interpreted 
as approximate ground state solutions of the GPE of the 
left and right wells. The equations of motion for the rel- 
ative population r\ — (N 2 — N%)/N and phase 8 = 6 2 — 9\ 
between the two symmetric traps are 



hr,= {2Ej/N)y/l^f sin {&) , 
2Ej r, 
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where E c — 2g J dr if>i (r) 4 is the variational ana- 
log of the capacitive energy in SJJ, while Ej = 
-N J drip i (r) [H q + gNi/jj (r)] tp 2 (r) is the Josephson 
coupling energy. The current-phase relation I = 
I c y/l — rf sin(0) is directly related to Eq. (||) where 
the critical current is given by I c — Ej/h. 
F l z (t) represents the contribution to the chemi- 
cal potential difference in the two wells due to 
the laser displacement l z (after linearizing in l z ), 
and where F = $ dr (ipi(v) 2 - ip 2 {v) 2 ) -§-Vi aser ~ 
Jdrz (V>i(r) 2 - V^r) 2 ). Theab ove variational 
method provides a simple and useful interpolating 
scheme between the low interacting limit N 2 E C <C Ej 
and the opposite limit N 2 E C 3> Ej. In the last case, and 
with t] <C 1, we recover the driven-pendulum phase equa- 
tion ||) and the critical velocity relations (||) and (JsJ) with 
Tiijp = y/Ej E c . In particular, it is legitimate to consider 
the Josephson coupling as a perturbation, with the the 
phase dynamics entirely determined by the difference in 
the chemical potentials /xi(iVi,Z z ) and /i 2 (N 2 ,l z ) in the 
two wells. In this case E c corresponds to 2 (dfix/dNi) l 
and h 2 uj 2 = Ej (d^i/ dNi) l . The critical velocity is pro- 

portional to the critical cu.ent: v c = (^) ' I m with 
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and {dn\/dl z ) N being F/2 in Eq.(^). These deriva- 
tives can be computed numerically. In the Thomas- Fermi 
(TF) limit they reduce to 
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where Vtf is the volume of the region in which is 
different from zero (in the TF approximation). 

We make the comparison of Eqs. (||) and (^) with a full 
numerical integration of the GPE in an experimentally 
realistic geometry relative to the limit N 2 E C ^> Ej. In 
particular, we show that Eq. (^), derived in the limit 
of rj -C 1, still remains a good approximation even for 
?y sa 0.4. The details of the numerical calculation are 
given elsewhere H|. 

We have considered the JILA setup, with N = 5 x 10 4 
Rb atoms in a cylindrically symmetric harmonic trap, 
having the longitudinal frequency loq — 50 s _1 and the 
radial frequency ui r = 17.68 s _1 . The value of the scatter- 
ing length considered is a — 58.19 A. A Gaussian shaped 
laser sheet is focused in the center of the trap, cutting it 
into two parts. We assume that the (longitudinal) 1/e 2 
half- width of the laser barrier is 3.5 /mi and the barrier 
height Vo/h = 650 s" 1 . 

Although the lifetime of a trapped condensate can be 
as long as minutes, we have made a quite conservative 
choice, by considering a time scale on the order of one sec- 
ond. The possibility to perform experiments on a longer 
time-scale will improve the observability of the phenom- 
ena we are discussing. With this choice of time scale, 
that corresponds only to few plasma oscillations, an adi- 
abatic increase of the velocity is not possible, therefore 
we proceed as follows. For t < the laser is at rest in the 
middle of the trap, l z — 0, and the two condensates are 
in equilibrium. For t > the laser moves across the trap, 
with constant velocity, and the relative atomic popula- 
tion is observed at i/ = 1 s. With this initial condition, 
which introduces small plasma oscillations in the rela- 
tive population, it is expected, in absence of dissipation, 
to slightly reduce the critical current by the numerical 
factor ss 0.725 (see the general properties of the driven 
pendulum equation Q). 

In Fig.l we show the relative condensate population 
V = (N2 — Ni)/N, calculated after 1 second, for differ- 
ent values of the laser velocity v. The crosses are the 
results obtained with the full numerical integration of 
the time-dependent GPE (Q) . The dot-dashed line shows 



the equilibrium values r/ eq of the relative population cal- 
culated with the stationary GPE and with the laser at 
rest in the "final" position l z = v tf. The displacement 
of rj (tf) from r/ eq is a measure of the chemical potential 
difference, being A^i — /i 2 — ^\ ~ NE c (r](tf) — r] eq )/2. 

For v < 0.42 /xm/s, the atoms tunnel through the bar- 
rier in order to keep the chemical potential difference 
A/i locked around zero. The dc component of the tun- 
neling current is accounted for by an averaged constant 
phase difference between the two condensates. This is the 
close analog of the dc Josephson effect in superconduct- 
ing Josephson junctions. The small deviations between 
the dashed line and the crosses are due to the presence 
of plasma oscillations (induced by our initial condition). 
At v 0.42 fim/s there is a sharp transition, connected 
with the crossover from the dc-branch to the ac-branch 
in SJJ. For v > 0.42 [im/s, the phase difference starts 
running and the population difference, after a transient 
time, remains on average fixed. A macroscopic chemical 
potential difference is established across the junction. In 
this regime ac oscillations in the population difference are 
observed. The frequency of such oscillations are approx- 
imatively given by A[i(t)/h (not visible in the figure). 
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FIG. 1. Fractional population imbalance versus the veloc- 
ity of the laser creating the weak link. A sharp transition 
between the "dc" and the "ac" branches occurs at a bar- 
rier critical velocity. The solid line and the crosses are the 
analytical and the numerical calculations, respectively. The 
dashed-dot line represents the static equilibrium value r] sq cal- 
culated with the center of the laser at vtf. 

The solid line of Fig.l corresponds to the solutions of 
Eqs. (H) and @ in which the value of the energy integrals 
E c N/h = 2.46 mr 1 and Ej/Nh = 2.41 x 10~ 4 ms" 1 
are chosen in order to give the correct value of lj p = 
2.44 x 10~ 2 ms- 1 and I c = 12.1 mr 1 . The val- 
ues uip, I c are calculated numerically studying the fre- 
quency of small oscillations around equilibrium and the 
current-phase relation, respectively. The force integral is 
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F/h — 1.060 ms" 1 \rni~ 1 . The parameters uj p , I c and F 
are calculated with the laser at rest (v — 0) in l z = 0. 
Using these values in Eq. (||) and taking into account the 
reducing factor 0.725 we obtain the value 0.407 \im s _1 
for the critical velocity, in agreement with the value ob- 
served in the simulation. 

Small deviations between the variational solutions (full 
line in Fig.l) and the numerical results (crosses in Fig.l), 
above the critical velocity, are due to "level-crossing" ef- 
fects. Numerical results show that when the con- 
densate ground state of the "upper" well is aligned with 
the excited collective dipole state in the "lower" well, 
a finite number of atoms go from the "upper" well to 
the "lower" well. Close to this tunneling resonance it is 
possible to control, by manipulating the barrier veloc- 
ity below a fraction of v c , the dc flux of atoms from the 
ground state condensate in the "upper" well to the longi- 
tudinal intrawell collective dipole mode of the condensate 
in the "lower" well. This effect is directly observable in 
the macroscopic longitudinal oscillations of the two con- 
densates (at frequencies ~ u>q). 

Concerning a possible realization of the phenomenon 
described in this work, we note that for small barrier 
velocities v, the motion of the laser sheet with respect 
to the magnetic trap with velocity v or, viceversa, the 
motion of the magnetic trap with velocity —v, are equiv- 
alent, there being negligible corrections due to different 
initial accelerations. 

Thus far we have discussed the zero temperature limit. 
At finite temperature dissipation can arise due to inco- 
herent exchange of thermal atoms between the two wells. 
This can be described phenomenologically by including a 
term —E c G9/uip in Eq. (^) where G is the conductance. 
Dissipation will be negligible as long as the characteristic 
time scale (EcG)^ 1 « (20G/K) s is bigger than the time 
scale of the experiment (~ Is). 

To conclude we note that while it could be difficult 
to measure directly the plasma oscillations, since their 

amplitude is limited by Arj < jq\J^f-, the macroscopic 
change in the population difference may be easily de- 
tected with standard techniques. Moreover the frame- 
work that we have discussed can be easily adapted to 
investigate the internal Josephson effect. 

Our phenomenological equations are similar to the 
driven pendulum equation governing the Josephson ef- 
fects in SJJs. As a consequence, within this framework 
we can study the "secondary quantum phenomena" , such 
as the Macroscopic Quantum Tunneling between differ- 
ent local minima of the washboard potential (see for in- 
stance fl9t|). 
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